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Cavitation is a limiting phenomenon in many domains of fluid mechanics. Instabihties 
of a partial cavity developed on an hydrofoil, a converging-diverging step or in an inter- 
blade channel in turbomachinery, have already been investigated and described in many 
previous works. The aim of this study is to evaluate a passive control method of the sheet 
cavity. According to operating conditions, cavitation can be described by two different 
regimes: an unstable regime with a cloud cavitation shedding and a stable regime with 
only a pulsating sheet cavity. Avoiding cloud cavitation can limit structure damages since 
a pulsating sheet cavity is less agressive. The surface condition of a converging-diverging 
step, like a Venturi-type obstacle, is here studied as a solution for a passive control of 
the cavitation. This study discusses the effect of an organized roughness, in the shape 
of longitudinal grooves, on the developed sheet cavity. Analyzes conducted with Laser 
Doppler Velocimetry, visualisations and pressure measurements show that the grooves 
geometry, and especially the groove depth, acts on the sheet cavity dynamics. Results 
show that modifying the surface condition, by varying the grooves geometry, can reduce 
cavity sheet length and even suppress the cloud cavitation shedding. 
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1. Introduction 

Cavitation is a crucial phenomenon encountered in fluid mechanics and for instance in 
turbomachinery domain. The unsteady character of the sheet cavity behavior on suction 
side of hydrofoils, on converging-diverging obstacles or on blades in turbines and pro- 
pellers is responsible for many issues like erosion, noise and vibrations. However, in many 
industrial devices, cavitation inception can not be avoided, so the challenge consists in 
a better understanding of the two-phase flow dynamics once cavitation has appeared, in 
order to reduce its negative effects. Simple geometries like 2D foil sections, or converging- 
diverging steps lik e Venturi-type sections are used to approach sheet cavitation dynamics 
( Barre et al. I (2009l) . A low pressure zone appears in the flow, downstream this type of ob- 



stacles. When this pressure becomes lower than the fluid saturation pressure, cavitation 

is induced and and a two-phase flow is obtained. 

Sh eet cavity dyn amics can be des cribed wit h a well- known cyc le ("KnapplllQSSl: lAvellan et al 



1991: Stutz fc Reboud 1997: Callenaere et al. 2001: Dular et al. 2004; Coutier-Delgosha et al\ 
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Sheet cavity grows from the Venturi throat until a re-entrant jet appears in 



the sheet cavity closure, and flows upstream, near the wall, below the cavity. The re- 
entrant jet separates then the sheet cavity when it reaches the liquid-gas interface of 
the flow. A large vapour cloud is shed and convected downstream the flow, while the 
sheet cavity length is substantially reduced. Then, the cloud of vapor collapses in a 
higher pressure region and sustains the re- entrant jet due to cloud implosion, which later 
produces a new shedding o f vapor cloud ( De Lange fc De Bruin 1997 : Callenaere et all 
200 ll: iBergerat et al. 2012). This cycle is characterized by the shedding frequency /. 
Such oscillatory behaviour of the sheet cavity is a witness of two-dimensional and three- 
dimensional cavity instabilities. These instabilities i nduce vibrations which are a disabil- 
ity for many industrial applic ations, as for inducers ( Bakir et al. 2003 : Mejri et al. 20061 : 
ICampos-Amezcua et aLll201o( ). 

Many studies have analyzed this sheet cavity cycle observed on cavitating hydrofoils 
or V enturi profiles. They have resulted in th e definition of two types of oscillatory behav- 
iors (|Callenaere et aLll200ll : [Savvaadill2010[ ): the first one called cloud cavitation regime 
appears when the sheet cavity is long enough and generates periodical shedding as it has 
been described previously. But when t he sheet cavity length is shorter because of a higher 
pressu re level (Callenaere a/.ll200ll) and/or a Reynolds number smaller than a critical 



value (|Keil adl201lU2012[) ." no cloud cavitation shedding occurs. In this regime, the 



sheet cavity pulses but the cavity closure is always at the same location and no large scale 
detached cavitation structure is observed. In order to characterize these two regimes, a 
Strouhal number StL is defined (as it is presented in section [231) ■ For a sheet cavity 
regime, StL —0.1 while a cloud cavitation regime is characterized by 0.2 < St^ < 0.4. 
The cloud cavitation regime is more aggressive and leads to increased damages of the 
solid structures.. It is thus interesting to find a way to control the sheet cavity instabil- 
ity in order to limit erosion and/or noise. A passive control method of the sheet cavity 
regime has therefore to be evaluated by its ability to reduce the sheet cavity length, in 
order to favorise nearly stable regimes in front of periodical unsteady behaviors. 

Some studies have already been conducted to delay cavitation inception by modifying 
the foil surface rough ness using local protuberances wi th different geometries or dis- 
tributed irregularities ( Arndt fc Ippen 1968 : ArndtlfigSl ). This modification acts on the 
turbulent boundary layer which drives the onset of sheet cavitation. However, in many 
industrial applications, the effect of developped cavitation on surface drag and/or per- 
formance is also of primary importance. Therefore, attention should be paid not only 
on delaying cavitation inception, but also on moderating fully developped cavitation. 
Only a few studies exist about the role of the surface condition in a fully developed 
cavitating flow. Some of these works show that roughness is able to decrease sheet cav- 
ity leng th, to increase oscillation frequency of the cavity or even to ch ange cavitation 
regime (jKawanami et aZ.lll997l : IStutzll2003t ICoutier-Delgosha et aZ.ll2005l ) . The roughness 
distribution and geometry is a major characteristic of this passive method of cavitation 
control: tranversal or longi tudinal grooves, wit h smooth or straight edges, for different 



depths or interval lengths ([Yoniian et al. 2010f ). 



The present study investigates the effect of different grooved suction side surfaces of 
a Venturi-type section with convergent and divergent angles respectively 18° and 8°. 
Different grooved plates have been studied in order to emphasize the crucial geometric 
parameter of this organized roughness to obtain a passive control of cavitation. High 
speed visualizations of the cavitating flow show the effect of flow modifications near the 
wall on the sheet cavity development, and next, velocity measurements in non-cavitating 
conditions clarify the effect of the grooved surface on the flow structure. 
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2. Experimental set-up 

Experiments were conducted in the closed loop test rig of the DynFluid laboratory wa- 
ter tunnel. This rig is composed of two storage tanks with a capacity of 4m'^ each. The 
centrifugal pump can reach 1450 rpm. The test-section for experiments is presented in fig- 
ure[T] The study flow volume, which is presented in figurefU measures 120 x 100 x 800 mm^ . 
The bottom wall is made of a convergent (18° angle) and a divergent (8° angle), which 
results in a minimum height Hthroat = 67 mm at the Venturi throat. The flow rate is fixed 
for all measurements at Q =240m^.h~^, which gives a maximun velocity at the Venturi 
throat Vthroat — 8m.s~^. Considering that the water viscosity in the operating condi- 
tions is v =1.007 10~^m^.s~^, then the Reynolds number is Rethroat= VthroatHthroat/v ~ 
5.5 10^, for the water at a temperature T = (19 ± 1)° C. Several honeycombs followed 
by a circular contraction provide the flow a velocity top hat profile upstream the Ven- 
turi with an inlet velocity equal to Vref =5.56 m.s~^, for all experimental configurations, 
with less than 3% of turbulence intensity (figure [2]), measured at a distance y =300 mm 
upstream the Venturi throat (at an undimensional distance y* — y / Hthroat — — 4.5 from 
the throat). A vacuum pump can decrease the pressure in the test rig by decreasing pres- 
sure at the free surface in each storage tanks in order to acts on the cavitation number 
c = {Pref — Pvap)/ i^P'Vref'^), with Pref the pressurc measured at y =250mm upstream 
the Venturi throat (at y* = y/ Hthroat — — 3.73). The test section includes four plexiglas 
windows to permit three-dimensional visualisations of the flow. 

The bottom wall of the test section, in the divergent part of the Venturi, is made of 
interchangeable plates of dimensions 7 mm thick, 242 mm long, and 120 mm wide, to cover 
the entire Venturi surface downstream the Venturi throat. The junction between the plate 
and the Venturi basis is done upstream the Venturi throat to minimize its disturbances 
on the flow. So disturbances are expected to be smoothed before the flow comes on the 
throat, thanks to the favorable pressure gradient. Seven 2 MI-PAA KELLER absolute 
pressure sensors even out with the Venturi basis wall and distributed on the Venturi 
basis surface (figure [3] and table [2|). These sensors acquire pressure measurements during 
1 min at a sampling rate of 1000 Hz. The plate is screwed on the basis by ensuring the 
sealing between both components with a vacuum grease. Table [1] presents the different 
grooved plates used for experiments. Plate will be considered as the reference case, 
with a smooth surface on the suction side of the Venturi. The other plates, from 1 to 8, 
have grooved surfaces. These grooves are made with a 1mm or 2 mm diameter bur d, 
with a h depth and the width of the ridge is e. The distance between the middle of two 

adjacent grooves is equal to A = e+2^/ dh — for plates 1 to 6 and A = e -|- d for plates 
7 and 8. Grooves open to the throat, so roughness effects already from the onset of the 
sheet cavity. All positions in the test section are expressed with undimensional values 
X* = x/Hthroat, V* = v/Hthroat and z* = z/Hthroat, and the origin of the cartesian 
coordinate system is located in the throat, in the middle of the test section width. 

The test of these grooved plates and their comparison with the results obtained with the 
smooth reference plate allow to evaluate the effects the effect of the surface condition of 
the Venturi bottom wall on the flow dynamics near the wall, with or without cavitation. 



3. Effect of grooved surfaces on cavitating flow 

Visualisations are conducted for each plate (the smooth and the grooved ones), for dif- 
ferent values of cavitation number a defined as cr = {pref — Pvap) / {\p^ref'^)- In this defi- 
nition, Pref and Vref rcprcscnt the inlet pressure and velocity of the test section, measured 
250 mm upstream the Venturi throat, while pvap is the vapor pressure at temperature 
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Figure 1. Test section of the experimental device using a Venturi profile (the origin of the 
cartesian coordinate system is located at the throat, in the middle of the test section wide, 
Q =8° is the angle for the LDV measurements system xy' z' and Hthroat =67 mm is the reference 
length) . 
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Figure 2. Characterisation of the flow in the test section inlet: (a) undimensional longitudinal 
velocity profile Vy = Vy/vref (the line represents the mean value of the velocity Vref =5.56 m.s^^) 



and (6) the profile of the turbulence intensity 



RMS(Vy) 



where < Vy > is the time-averaged 



velocity, y* 



4.5 (Results come from LDV measurements, presented in the section UjT} . 
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Figure 3. Venturi components: (a) Venturi basis with CI to C7 pressure sensors on the basis 
surface, (6) junction between the Venturi basis and a plate and (c) zoom of a grooved plate 
with the definition of the grooves geometric parameters {d and h are respectively the grooves 
diameter and depth, e is the width of the ridge and A is the grooves wavelength). 



T —292 K of the fluid with a density p. Temperature deviation during measurements is 
small enough (about 1 K) to neglect temperature effect on results but in order to have 
precise values of cr, the flow temperature is measured at each acquisition (with an un- 
certainty of 0.1 K). Temperature for all measurements is then considered as the ambient 
temperature and pvap — 2200 Pa. For each studied configuration, a ranges from 0.9 to 
1.7 with 1% of uncertainty. 
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Table 1 . Characteristics of the studied plates used on the Venturi suction side 
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Table 2. Pressure sensors on the Venturi suction side 



3.1. Sheet cavity visualisations 

Visualisations of the sheet cavity developed on the Venturi consist in the acquisition of 
4000 images for each configuration (for a given cavitation number and a given grooved 
plate geometry). A CamRecord 600 Optronics camera is used with a Zeiss Makro-planar 
T* which has a focal length of 100 mm. This high-speed CMOS camera records images 
with a 1280 x 512 pixels resolution, a rate of 1000 fps and an exposure time equal to 
0.167 IIS. The pixels size is 12 fim x 12 /im for an active area of 15.35 mm x 12.29 mm. The 
flow is illuminated by a Super Long Life Ultra Bright (SLLUB) White Led Backlight from 
Phlox on a 200 mm x 200 mm light output area. Its minimal luminance is 3000 cd.m~^ 
in continuous mode and its uniformity is 99.54%. 

3.2. Effects of the Venturi surface condition on the sheet cavity size 

In order to determine the cavity mean length, mean values and standard deviation of 
the grey levels are calculated with the 4000 images recorded in each flow configuration, 
in the (yOz) plane (side views of the attached sheet cavity dynamics). The instantaneous 
images are normalised by a reference image. The mean cavity length extends from the 
Venturi throat to the cloud cavitation separation location. This closure of the sheet 
cavity is detected by lo calising the maxim um value of the standard deviation of the 
grey levels according to Dular et al. ( 2004[ ) (figure |3|). Figure [5] shows the evolution of 



the cavity mean length L according to the cavitation number cr, for all plates. The 
uncertainty of L is less than 2% (by evaluating the effect of images processing) and a is 
calculated with 1% of uncertainty. Plates 3, 4 and 5 provide a sheet cavity length close 
to the the one obtained with the reference plate 0. Conversely, plates 0, 2 and 6 lead to 
significant differences of the cavity lengths. The results obtained with these three plates 
are thus detailed hereafter. It can be observed that for the plate 2, the cavity mean length 
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decreases as it has been enco untered in litterature, for random distributed roughness 
(|Coutier-Delgosha et a^.l 120051 ). This plate has a large number of grooves {N =124 and 
A =0.97 mm) and is characterized by a small depth h =0.25 mm. But other types of 
grooves have the opposite effect, i.e. an increase of the mean length of the sheet cavity, as 
we can see for the plate 6, which has a small number of grooves (A'^ =57 and A =2.1 mm) 
with a large depth /i = 1 mm. This difference between the effects of surface condition 
is larger for small cavitation numbers. When a is greater than 1.2 i.e. for small sheet 
cavities, the influence of the surface condition on the mean cavity length seems to be less 
important. The major difference between plates 2 and 6, is the depth parameter h, the 
diameter d and the number of grooves on the bottom wall of the Venturi. In comparison 
with plates 3, 4 and 5 which do not affect significantly the sheet cavity length with 
respect to the reference plate 0, the depth h seems to be the criterium that drives the 
change of effects on the sheet cavity. Indeed, the plate 2 is one of the plates which has 
the smallest depth. Plate 1 has an even smaller depth, but this plate seems to have no 
effect on the cavitation dynamics. The depth of plate 1 has been chosen in order to 
be of the order of the viscous sublayer thickness 6y (zone of the boundary layer where 
the velocity increases linearly from the wall, beyond this zone, the velocity increases 
logarithmic ally with the height). If the bottom wall of the Venturi is considered as a 
plane plate (jSchlichting fc Gerstenll2004l) . Sy is defined by: 



where Rcy* is the critical Reynolds number which determines the distance y* from the 
Venturi throat where the boundary layer becomes turbulent {Rcy- = 5 10'"' for a plane 
plate) and (7/ is the coefficient of friction expressed as: 



Ct 



K 



logi?ej 



-G {log Rey.) 



(3.2) 



with K = 0.41 the Karman constant and G a function that gives G (log [Rey)) = 1.5 
in the studied zone 10^ < Rcy, < 10^ (|Schlichting fc Gersten|[200l . 

In our experimental set-up, Sy is inferior to 0.20 mm {Sy / Hthroat — 0.003) at the end 
of the plate (at y* ~ 3.4). 



Indeed, results obtained with plate 1 (with a grooves depth inferior to the viscous 
sublayer thickness 6y) are identical to the one measured with plate 0. Conversely, if the 
depth of the grooves is too large, then the sheet cavity length increases. 

However, if the depth is increased up to /i = 2 mm, results are different: for plates 7 
and 8, we can see a decrease of the sheet cavity length, like for the plate 2 which presents 
a small depth. These two plates 7 and 8 are special because the geometry of the grooves 
is different. Grooves for plates 1 to 6 are arcuate hollows while grooves are cylindrical 
gutters for plates 7 and 8 (ft > This geometry leads to a decrease of the sheet cavity 
length. This phenomenon may be related to the fluid flow inside the grooves, which may 
be different with plate 6, compared with plates 7 and 8. Indeed, the number of grooves 
N is the same for these three plates. It can thus be deduced that this parameter is not 
crucial for the effect of grooves on the sheet cavity length. Only the depth h is different 
from plate 6 to plates 7 and 8. It implies that a reduction or an increase of the sheet 
cavity length can both be obtained with appropriate grooves, depending on their depth h. 
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Figure 4. Sheet cavity on the smooth plate for a =1.18: (a) normahzed average image, (6) 
root mean square of normahzed images (the line represents the undimensional sheet cavity length 
L* = L/Hthroat — 1.62) and (c) a profile of the root mean square of normalized images RMS(/) 
plotted in z*=0. 




Figure 5. Undimensional sheet cavity length L* = L/Hthroat according to the cavitation number 
a: X: plate 0, A: plate 1, <(>: plate 2, -A': plate 3 , +: plate 4, plate 5, •: plate 6, plate 7, 
□ : plate 8. 



The former analysis was focused on the effects of the grooves on the sheet cavity 
length. In addition, the effects on the mean height of the sheet cavity are also in- 
vestigated hereafter: for that purpose, the maximum height of the sheet cavity is de- 
tected. The calculation of the mean value of all these maximum heights (which can 
variate in position) gives information about the shape of the sheet cavity,, according 
to the grooved surfaces. Figure [6] displays the undimensional sheet cavity mean height 
H* = H / Hthroat (with 2% of uncertainty) according to the undimensional sheet cavity 
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Figure 6. Undimensional sheet cavity mean height H* = H / Hthroat according to the undimen- 
sional sheet cavity mean length L* — L/ Hthroat- X' plate 0, A: plate 1, <C>: plate 2, *: plate 3 , 
+ : plate 4, plate 5, •: plate 6, *: plate 7, □: plate 8. 



mean length L* = L/ Hthroat, which depends on the cavitation number, for all different 
studied plates. It can be observed that results obtained with plate 6 are significantly 
different from all other data. As grooved plates introduce a three-dimensional geometry 
with different grooves wavelengths, these Venturi surfaces induce a changing in the sheet 
cavity dynamics: the grooves wavelength forces the flow, and especially the inception of 
the sheet cavity, in the grooves hollow, but also the development of the sheet cavity or the 
re-entrant jet. This phenomen, with a three-dimensional wavelength imposed to the flow 
is diffe rent from the sm ooth case, with the plate 0. The cavity closure line is actually in- 
clined ( De Lang3 1996) and the cloud cavitation shedding is linked to three-dimensiona l 
instabilities ( Laberteaux &: Ceccio 1998 : Duttweiler &: Brennen 1998 : Foeth et aLlbOOSi ). 
three-dimensional effects are then to be taken into account to evaluate the efficiency of 
the passive control. As each plate presents a different spanwise distribution, we can say 
that depth is a crucial parameter but that the grooves wavelength A is also important. 
Plates 6, 7 and 8 have largest wavelength (A > 2 mm and N < 60). For small sheet 
cavities characterised by < L* < 1, which corresponds to 1.2 < cr < 1.9, all curves are 
nearly superimposed, which suggests that grooves play a less crucial role. 

Figure [7] which presents the aspect ratio H / L according to the cavitation number for all 
grooved plates, shows that most of grooved plates present large variations of the aspect 
ratio values with different cavitation numbers while the reference plate case leads to 
a nearly constant aspect ratio. The grooved plates 3, 4 and 5 for which the sheet cavity 
length is almost equal to the reference case, for all cavitation numbers, are characterized 
by an aspect ratio 0-1 < < 0.2. These values are still close to the reference case. On 
the other hand, the plate 1, with a grooves depth h too small to modify the sheet cavity 
length, has an aspect ratio very different from the smooth plate 0. This result shows that 
the grooves have effects not only on the sheet cavity length but also on its shape. Plates 
7 and 8 present an evolution of the aspect ratio H/L similar to the plate 1. The most 
obvious observation is that the evolution is inversed for plate 6. Indeed, for plate 6, the 
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Figure 7. Aspect ratio H/L according to the cavitation number a: x: plate 0, A: plate 1, <0: 
plate 2, *: plate 3 , +: plate 4, plate 5, •: plate 6, *: plate 7, □: plate 8. 



ratio H/L is high for 1 < tr < 1.4 and smaU for 1.4 < tr < 1.9 (0.2 < f < 0.4 for the 
first range of a and 0.1 < ^ < 0.2 for the second), while, for plates 7 and 8, H/L is 
first in the range of 0.1 < f < 0.2 when 1 < tr < 1.4 and then 0.2 < f < 0.4 when 
1.4 < a < 1.9. Cavitation dynamics is thus completely different for plates 6 and other 
plates, which have a smaller grooves wavelength A, excluding plates 7 and 8 for which 
the grooves geometry is quite different. It confirms that the crucial parameter in the 
action of the grooves is i) the depth, which may modify the re-entrant jet progression, 
ii) the wavelength, which can force the cavitation dynamics in the spanwise component. 
As it was remind ed in the introduction , the cavity regime depends on the mean size of 
the sheet cavity (jCallenaere et aLll200lh . An analysis of the Strouhal number based on 
the characteristic frequency of the sheet cavity is conducted to determine the effects of 
grooved surfaces on the unsteady behavior of the cavity. 



3.3. Two cavitation regimes 

The fluctuations of the two-phase area in the sheet cavity regime or the shedding of the 
cloud of vapor in the cloud cavitation regime are both periodical. The frequency of these 
phenomena is determined by studying the evolution of the grey levels at the sheet cavity 
closure position. The power spectral density of this sheet cavity length evolution allows 
to determine the characteristic frequency of the sheet cavity according to the mean cavity 
length for each studied case of Venturi profile grooved surface (figure HI), with 0.2% of 
uncertainty. 

Here, the frequency difference between the smooth Venturi surface and the other 
grooved plates is greater when the mean cavity length is large (it means for small cav- 
itation numbers). It can be observed in the range < L* < 1, which corresponds to 
1.3 < CT < 1.9, two peaks on the frequency spectrum, with similar amplitudes, at 43.95 Hz 
and 8.79 Hz. These two frequencies are harmonics and are both present in this studied 
range, as it is shown in figure [9] for the plate 0. This frequency of 43.95 Hz is the blade 
passing frequency of the centrifugal pump used to impose the flow rate: as the pump 
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Figure 8. Frequency of the sheet cavity length osciUation fi according to the undimensional 
sheet cavity mean length L*: x: plate 0, A: plate 1, <D: plate 2, plate 3 , +: plate 4, plate 
5, •: plate 6, *: plate 7, □: plate 8. 




Figure 9. Frequency spectrum of the variation of the gray level in the closure of the sheet 
cavity, on the plate 0: (a) for a =1.17, (b) for a =1.44 



with 5 blades operates at 528 rpm, the blade passing frequency is 44 Hz. It means that 
smallest cavities oscillate then at the blade passing frequency. For largest sheet cavities, 
so for a cavitation number a in the range 1 < it < 1.3, there is only one peak on the 
frequencies spectrum (figure IHl (a) presents results for a =1.17). This peak corresponds 
to the shedding frequency of the cloud cavitation and it decreases when the sheet cavity 
length increases. Plates 7 and 8 lead to a different feature for the frequency derived from 
the grey level analysis in the closure of the sheet cavity. The cavitation number a has no 
effect on the frequency peak of spectra, which is always /l =43.95 Hz. This result shows 
again that the cavitation dynamics is particular for these two plates 7 and 8. Figure [TOl 
shows frequency spectra obtained for plate 7 (figure Unta)) and plate 8 (figure rTOT b')') for 
a =1.17. 
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The characteristic frequency of the sheet cavity is also represented on the figure [TT] 
with the undimensional Strouhal number St^ = Lfi^/vthroat, with 6% of uncertainty. If 
we look at all plates excluding plates 7 and 8, we can observe that for 1 < a < 1.3, 
th e Strouhal number is a. r ound 0.3, as it has been reported in previous experiments 
in Coutier-Delgosha et al. ( 20061 ). This range of cavitation numbers leads to a cloud 



cavitation regime, characterized by large fluctuations of the cavity closure position and 
large cloud cavitation shedding. On the other hand, for 1.3 < cr < 1.9, the Strouhal 
number decreases until reaching a value smaller than 0.15. This range corresponds to a 
sheet cavity regime, for which there is no convected cloud cavitation, the cavity closure 
position is almost at a constant distance from the Venturi throat. In the sheet cavity 
regime, the sheet cavity pulses. These two types of variation of the Strouhal number 
exist for grooved p lates but also for the smooth ori e (plate 0), and have been brought out 
by other authors ( Far hat 1994 : Dular et allhOl^ . It can be also observed in figure [Til 
that the grooved plates surfaces can modify the Strouhal number. Indeed, smooth surface 
plate presents Strouhal numbers larger than grooved plates. 

Plates 7 and 8 have a quite constant Strouhal number 0.05 < Sti < 0.15, which is 
characteristic of the sheet cavity regime, with a pulsating sheet cavity and no cloud 
cavitation shedding. Figure [T^] compares the sheet cavity evolutions at a same tr =1.17, 
for the reference smooth plate and the grooved plate 8. It can be seen that in the case 
of the plate 8 the sheet cavity pulses (figure fT2la)) without any shedding contrary to a 
smooth surface of the Venturi profile (figure [T2lb)). Then, the grooves depth h controls 
the cavitation regime and can suppress the cloud cavitation for a large range of a. Plates 
7 and 8 are thus able to stabilize the sheet cavity and to reduce its length. It is thus 
possible to obtain a passive control by using specific geometrical parameters of grooves. 

With the study of the sheet cavity height evolution (section r3.2|) . we can also determine 
the frequency /h of the sheet cavity height. Figure [T3] compares this frequency /h with 
the frequency of the sheet cavity length evolution for plates to 6. Results for plates 
7 and 8 are not presented here because for all the studied range of a, there is no cloud 
cavitation shedding. For all grooved plates and for the smooth plate 0, the curve is linear 
so the frequency of the height evolution is linked to the length evolution frequency. 
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Figure 11. Strouhal number StL = accordine; to the cavitation number a: x: plate 0, A: 

^throat ^ ' 

plate 1, <C>: plate 2, *: plate 3, +: plate 4, A: plate 5, •: plate 6, plate 7, □: plate 8. 



3.4. Pressure measurements 

In order to analyze the different cavitation dynamics detected with the grooved surfaces, 
pressure measurements have been acquired at the bottom wall of the Venturi. Relative 
pressure measurements have been first conducted with no velocity and no use of the 
vacuum pump to define the zero setting state for each sensor. All pressure values for 
the cavitating fiow are thus negative. On figure [T31 values of the time-averaged relative 
pressure are presented from the Venturi throat down to a distance j/'* =1.25, with dif- 
ferent cavitation number a, for plates 0, 2, 6, 7 and 8. This figure diplays the differences 
between sensors C4 and C5 (both placed at ?/'* =1 from the Venturi throat. Indeed, 
plates 6, 7 and 8 present large differences between these two sensors, while they are posi- 
tioned at the same distance j/'*, but at different positions in the span x*: C4 is situated 
in X* =0, in the middle of the test section, while C5 is located at station x* =0.4478, 
near the wall. This result shows that grooves of plates 6, 7 and 8 play a major role on 
three-dimensional sheet cavity dynamics, only in the range 1 < cr < 1.4 (figure [Til (a), 
(b) and (c)). This range corresponds to the cloud cavitation regime. For the sheet cavity 
regime (figure fTiT d)). this difference between sensors C4 and C5 disapears. 

The analysis of the pressure fluctuations obtained with plates 0, 2, 6, 7 and 8, on the 
figure 1151 indicates that the maximum amplitud e of fluct uations of each pressure sensor 



corresponds to the sheet cavity closure dCallcna ere et aZ.ll2001il . For example, for a =1.5, 



the undimensional sheet cavity length is L* < 0.5 on the plate (figure [5]). The sensor 
CI is the only sensor located in the sheet cavity. Fluctuations of pressure measurements 
is thus an indicator of sheet cavity length. This figure confirms the results obtained with 
the image processing: if we look at the sensor CI, the maximum of pressure fluctuations 
is at a =1.5 for the smooth plate 0, a < 1.5 for the plate 2 and a > 1.5 for the plate 6, 
with larger grooves. This result shows that the length of sheet cavities is smaller on plate 
2 than on plate and a longer sheet cavity is obtained on plate 6 than on the others. As 
it was observed for the plate 8 previously, pressure fluctuations are a witness of the cavity 
regime changing. Indeed, in the figure ITST e). which represents results for the plate 8, no 




Figure 12. Visualisation of (a-j) the cloud cavitation shedding for the smooth plate and 
(a'-j') the sheet cavity pulsation for the plate 8, when a =1.17 (At = 2ms between two images) 
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Figure 13. Frequency of sheet cavity mean height variation fn according to the frequency 
calculated with the cavity length variation /l: x: plate 0, A: plate 1, <C>: plate 2, plate 3, +: 
plate 4, plate 5, •: plate 6. 



^ -6 




Figure 14. Time-averaged value of the relative pressue P according to the undimensional 
distance y'* — y/ Hthroat from the Venturi throat for (a) cr =1.18, (6) o =1.22, (c) cr =1.65 and 
(d) a =1.80: x: plate 0, 0: plate 2, •: plate 6, plate 7, □: plate 8. 



clear maximum of ^p-^ can be identified, so we are not able to detect the sheet cavity 
length by analyzing pressure measurements. This is an other indicator of the sheet cavity 
regime. The difference between a sheet cavity regime and a cloud cavitation regime can 
thus be demonstrated also by the analysis of the pressure measurements. 
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Figure 15. Pressure fluctuations RMS(P) on the Venturi bottom wall according to the cavita- 
tion number cr, at different distances from the ventrui throat:(o) for plate 0,(fo) for plate 2, (c) 
for plate 6, (d) for plate 7 and (e) for plate 8. A: CI, •: C2, □: C4, x: C6. 



4. Grooved surfaces effect on a non cavitating flow 

To discuss the reasons for the modifications of cavitating flows with the grooved sur- 
faces, the velocity profiles are measured in the non cavitating flow, and more especially 
in the boundary layer, as it is the cavitation inception zone. Therefore, velocity measure- 
ments with Laser Doppler Velocimetry have been conducted. 

4.1. Laser Doppler Velocimetry measurements (LDV) 

A one dimensional FlowExplorer Mini LDA system of Dantec Dynamics has been used 
to measure the longitudinal componant Vy of the flow. This system consists in a factory- 
aligned and calibrated optical probe, with a focal lens of 300 mm, and a signal processor. 
The optical head provides two laser beams with 25 mW power and a wave length 660 nm. 
The measurement volume dimensions are 0.1mm in diameter and 1mm in length. The 
measurable velocity fluctuation is 0.002% of the velocity range. This sytem has a high 
accuracy, as the calibration coefhcient uncertainty is lower than 0.1%. In each measure- 
ment position, 10000 samples arc acquired, with a limiting time of acquisition fixed of 
60 s. The data rate is near 300 Hz, while the validation rate is upper than 60%. For mea- 
surements downstream the Venturi throat, the optical head is inclined at 8° in order to 
align the laser beams with the Venturi slope. In this way, measurements can be acquired 
very near the wall so velocity measurements can be used to analyze the effect of surface 
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condition on the boundary layer of the flow. Then, the acquired velocity component is 
not the longitudinal velocity of the test section but the velocity component parrallel to 
the slope of the Vcnturi bottom wall. In the discussion of LDV results, this component 
is then called Vy in order to clarify explanations. 

On the other hand, velocity measurements presented on figure [2] were made with a 
horizontal optical head, so the component is really the longitudinal velocity Vy in this 
case. 

4.2. Effects of the grooves on the velocity profile of the non cavitating flow 

Figure [T6l shows undimensional velocity profiles Vy* = v'y/v'^, where Vy is the flow ve- 
locity component parallel to the bottom wall Venturi slope and v'^ is the same velocity 
component measured far from the plate wall, at 2 * = z' / Hthroat =0.9 (this zone is not 
shown on the figure) . Velocity profiles are represented here at different distances from the 
throat, in y* =0.5, y* =1 and y* =2. We can see on this figure [T51 that grooves modify 
the flow, near the wall, up to z'* =0.07 for plates 2, 6 and 7 and until z'* =0.3 for plate 8. 
This modification depends on the geometry of the grooves. For all values of y*, excluding 
the plate 8 for which the velocity profile is totally different, the velocity is always lower 
for the plate 6, which is the plate with the larger grooves and minimum number N of 
grooves. On the other hand, for smaller grooves, with the plate 2 for example {d = 1 mm, 
h = 0.25 mm and e = 0.1 mm) which has a large grooves number N =124, in y* =0.5, the 
velocity profile is similar to the velocity profile of the reference plate 0, without grooves. 
At station = 1 or y* = 2 the velocity for the plate 2 becomes lower than for the 
reference plate 0. The effect of plate 2 on the grooves overlaps even the curve obtained 
with plate 6. Then the effect of the grooves geometry disapears from y* = 2, but grooves 
still play the role of a brake for the flow. Very close to the wall, from z'* = until 
z'* = 0.02, the velocity profile Vy* for the plate 7 is superimposed on the velocity profile 
of the plate 0, when < y'* < 2. This feature explains why plate 7 reduces sheet cavity 
length like plate 2. If z'* > 0.02, then Vy* is lower than other plates 2 and 6. On the other 
hand, the velocity profile obtained for the plate 8 is totally different. Vy* is very low from 
z'* = 0.3 for 0.5 < y* < 2. This analysis permits to bring out the effect of structured 
roughness surface on the flow but also the importance of the geometric parameters of 
the grooves. The grooves wavelength is thus a crucial parameter to infiuence the sheet 
cavity dynamics. 

In order to investigate in more details the effects of each groove on the velocity profile, 
two specific velocity profiles have been acquired for each plate: the first one is located in 
the middle of the test section wide, at x* = 0, above the recess of the groove, while the 
second one is measured at x* = A/2, above the edge of the groove. It can be observed 
(figure [T7| that the relative difference between the two curves is always lower than 3%. 
The effect of the grooves is thus not local but global, by the way it changes the fiow 
instabilities. 

This non-cavitating fiow analysis provides a basis for analysis of measurements per- 
formed in cavitating conditions. As it was reported previously, the use of plate 6 results in 
a decrease of the velocity fiow near the wall. The fluid encountered a larger surface to flow 
with larger grooves, so the velocity is considerably decreased. But with the analysis of 
the mean sheet cavity length, we can say that larger grooves increases the cavity length. 
Moreover, with larger grooves, the re-entrant jet can travel further upstream before reach- 
ing the interface and cutting the sheet cavity, in order to extract a cloud cavitation. The 
immediate consequence is that the mean sheet cavity length should be smaller with large 
grooves. An hypothesis that can be proposed is that some grooves provide small recircu- 
lation zones that prevent the re-entrant jet to flow upstream, underneath the sheet cavity, 
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Figure 16. Longitudinal non-dimensional velocity profile Vy * /v'y^ of the non-cavitant flow, at 
(a) y* — 0, (&) y* = 0.5, (c) y* — 1 and (d) y* = 2. x: plate 0, <C>: plate 2, •: plate 6, *: plate 
7, □: plate 8. 




Figure 17. Longitudinal non-dimensional velocity profile Vy * /vy^ of the non-cavitating flow, 
at y* = 0.5, for (a) the grooved sheet 6 and (6) the grooved sheet 2. •; x* = 0, A: x* — ^ 



as it is discussed in lCoustolsl (j200l[ ). for the study of a non-cavitant flow boundary layer. 
Then the re-entrant jet is stopped by these recirculation zon es. Not only re-entrant jet i s 
crucial for the cavitation shedding but also side-entrant jets ( Koop fc Hoeiiinakerj2010l ). 
The three-dimensional distribution of grooves and their geometry have of course effects 
on the side-entrant jets dynamics. 



5. Conclusions 

Effects of the surface condition of a Vcnturi profile have been investigated using visu- 
alisation, Laser Doppler Velocimetry, and pressure measurements. The aim of this study 
was to evaluate the importance of grooved surfaces provided by machining. Geometry 
of these grooved surfaces have been observed in order to identify crucial parameters. 
Results show that the depth of grooves is a determining factor. It is demonstrated that a 
grooves depth smaller than the viscous sublayer thickness has no effect on the sheet cav- 
ity length. Some plates lead to a cavitation instabilities changing: essentially plates 7 and 
8, for which the sheet cavity length is reduced and the plate 6 which increases this sheet 
cavity length. These three plates have a large depth /i > 1. If the depth is h > ^, then 
the grooved plate is able to reduce the sheet cavity length. The study of the shedding 
frequency of cloud cavitation, by image processing highlighted two cavitation regimes: 
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for small cavitation numbers (1 < ct < 1.3), cavitation is in an unstable regime, with 
large sheet cavity and large shedding of cloud cavitation. This regime is characterized 
by a Strouhal number St^ ~ 0.3. The second cavitation regime is a stable regime, with 
1.3 < cr < 1.9, which presents sheet cavities with only oscillations in the downstream 
part (this oscillation being related to the blade passing frequency of the circulation pump 
used in the test rig). A large grooves depth h, can modify the sheet cavity regime. Plates 
7 and 8, with h = 2 mm, are examples of roughness surfaces that can suppress the cloud 
cavitation shedding in a large range of cavitation number a. This study has thus demon- 
strated the feasability of a passive control of cavitation on a Venturi profile by modifying 
the surface condition of the bottom wall using distributed organized roughness. 
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